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Abstract: Consistent factorization theorems in high-energy scattering near the endpoint 
are presented in the framework of the soft-collinear effective theory. Conventional factor- 
ization theorem separates the soft and collinear parts successfully, but each part encounters 
infrared divergence and mixed ultraviolet and infrared divergences. We present factoriza- 
tion theorems in which the infrared divergences appear only in the parton distribution 
functions and the mixed divergence is removed by carefully separating and reorganizing 
collinear and soft parts. The underlying physical idea is to isolate and remove the soft 
contributions systematically from the collinear part in loop corrections order by order. Af- 
ter this procedure, each factorized term in the scattering cross sections is free of infrared 
divergence, and can be safely computed using perturbation theory. This factorization pro- 
cedure can be applied to various high-energy scattering processes. We show factorization 
theorems in Drell-Yan processes, deep inelastic scattering and Higgs production near the 
endpoint. 
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1 Introduction 

Theoretical predictions on high-energy scattering are based on the factorization theorem, in 
which the cross section is factorized into the hard, collinear and soft parts. It means that, 
though the strong interaction affects the scattering processes at different scales in various 
stages, these effects can be separated according to their kinematic regimes. The hard part 
consists of the partonic cross sections with hard momenta. The collinear part describes the 
effects of the collimated energetic particles. For hadron colliders, the collinear part involves 
the parton distribution functions (PDF) in the initial state. If an exclusive process is 
considered, say, a pion production, there appears the fragmentation function which depicts 
the hadronization process from a parton into a hadron. This is also classified as the collinear 
part. The soft part describes soft gluon exchanges among different collinear sectors. These 
three parts are decoupled and there is no communication among them. Therefore it is 
possible to compute the contributions for each part separately using perturbative quantum 
chromodynamics (QCD). 

The factorization property in various high-energy scattering processes has been well 
established in full QCD [1, 2]. In full QCD, once the scattering amplitudes are written, 
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they are examined closely by dividing collinear, soft momentum regions. A laborious 
analysis enables the factorization proof that collinear and soft parts are decoupled. The 
factorization proof becomes more elaborate in soft-collinear effective theory (SCET) [3-5] 
since the separation of collinear and soft parts is established from the outset at the operator 
level. The scattering cross sections are written as a product of the hard coefficient in terms 
of the Wilson coefficients, and the collinear and soft parts, which are expressed as the 
matrix elements of gauge-invariant operators. The factorization of the hard, collinear and 
soft parts in the hard scattering is now on a firm basis. We call this development the 
conventional factorization. 

Since the advent of Large Hadron Collider (LHC), precise theoretical predictions on 
high-energy scattering can be under experimental scrutiny. As well as the fixed higher-order 
corrections in the strong coupling a s , the resummed results using the renormalization group 
technique are available when there are large logarithms [1, 2, 6-9]. The theoretical accuracy 
available for the comparison with experiment can be attained by computing the factorized 
parts at higher orders or resumming large logarithms. The conventional factorization 
theorems have been successful because the collinear and soft parts can be expressed in 
terms of gauge-invariant operators and there is no communication between them. However, 
the conventional factorization theorems are plagued by infrared divergences. 

In massless gauge theories like QCD, there exist infrared (IR) divergence as well as 
ultraviolet (UV) divergence. When a massless quark emits a gluon, collinear divergence 
appears when these two particles are collinear. Soft divergence emerges when the energy 
of the emitted gluon becomes soft. If the emitted gluon is both soft and collinear to 
the quark, soft and collinear divergences occur simultaneously. This type of divergence 
is referred to as the IR divergence. The UV divergence also shows up, but they can be 
removed by counter terms. The remnant after removing the UV divergence governs the 
scaling behavior. On the other hand, the IR divergences cannot be removed by hand. 
Instead, physical quantities should be free of IR divergence. The Kinoshita-Lee-Nauenberg 
(KLN) theorem [10, 11] succinctly states that the IR divergences appearing in real gluon 
emissions and in virtual corrections at any given order in a s should cancel in scattering 
cross sections. 

There is no IR divergence in the hard part, and most of the IR divergences are cancelled 
in the sum of the soft and collinear parts. The exception is the radiative correction to the 
PDF. But since the PDF is basically a nonperturbative quantity, the IR divergence can 
be absorbed in the PDF. However, as we will explain, the soft and collinear parts in 
the conventional factorization theorems contain IR divergences, and furthermore mixing 
between IR and UV divergences. Therefore the collinear and soft parts themselves cannot 
be regarded as physical quantities. The evolution of the collinear and soft parts via the 
anomalous dimension does not make sense either because these parts contain IR divergences 
and mixed divergences, hence the renormalization scale dependence does not stem from 
UV divergences. 

Note that the IR divergences are not artifacts in calculating radiative corrections, and 
they always appear irrespective of the regularization methods. If we use the dimensional 
regularization to regulate both UV and IR divergences, the UV (IR) divergences appear as 
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poles of euv (^ir)- If the UV divergence is regulated by the dimensional regularization, and 
the IR divergence is regulated by the offshellness of external particles, the UV divergence 
still appears as poles in euv ; but the IR divergence appears as the logarithms of the 
offshellness. Unless the IR divergences are taken care of, the physical interpretation of 
each part with IR divergences is not possible, and we need more than the conventional 
factorization method. 

Our new idea of factorization respects the philosophy of the conventional factorization 
theorems in the sense that the cross sections are factorized and the collinear and soft 
parts can be expressed in terms of the gauge-invariant operators. But it goes further 
by reorganizing the collinear and soft parts in such a way that the IR divergence resides 
only in the PDF, as in full QCD, and the remaining parts in the factorization formula 
are free of IR divergences. Only after this reorganization, the IR-finite quantities can be 
computed systematically in perturbation theory. However, this new factorization is not 
just reorganizing the factorized parts to shuffle IR divergences as a matter of convenience, 
but there is a deeper physical reason. In loop calculations, the collinear part unavoidably 
covers the kinematic region with soft momentum when the loop momentum becomes soft. 
But the soft region is already accounted for in the soft part. Therefore the soft contribution 
in the collinear part should be removed in an appropriate way to avoid double counting. 

This removal can be systematically performed in SCET. The systematic technique 
to remove the soft contribution consistently from the collinear part is called the zero-bin 
subtraction [12] in SCET. The basic idea of the zero-bin subtraction is that, for each 
collinear Feynman diagram, the soft contribution can be computed by taking the limit 
where the loop momentum becomes soft, and then this soft contribution is subtracted 
from the naive collinear loop calculation. The zero-bin subtraction is not an esoteric 
technique used in SCET, but a simple tool to avoid double counting. The same idea can 
also be realized in full QCD. In full QCD, a careful separation of kinematic regions with 
different momenta should be performed, but the extraction of the soft contribution from 
the collinear part is tedious. SCET is an effective theory best fit to this purpose. SCET 
involves all the relevant fields necessary in the factorization, hence the factorization proof 
in SCET is equivalent to that in full QCD itself, but with more transparency. 

We emphasize that care should be taken to perform an appropriate zero-bin subtrac- 
tion. Naively any soft limit of the collinear loop momentum may do for calculating the 
zero-bin contributions, but in fact the soft limit of the collinear momentum should be de- 
termined by the size of the momentum entering the soft part. This is the essential point in 
introducing the zero-bin subtraction. If Q is the large scale in the system, and A is a small 
parameter in SCET, the soft part describes the fluctuations of order QX or QX 2 , and the 
soft limit in the collinear part should also be of the same order. If the soft part describes 
the interaction with momentum of order QX (QX 2 ), the soft limit in the collinear part 
should also be of order QX (QX 2 ). If the soft momentum in the soft part is of order QX, 
while the soft limit in the collinear part is of order QX 2 , or vice versa, there is a mismatch 
between the soft part and the soft contribution in the collinear part. With this mismatch, 
the zero-bin subtraction does not produce any sensible results, and the double counting is 
not removed completely. 



- 3 - 



After the double counting problem is handled appropriately, the scattering cross sec- 
tions near the endpoint again factorize, but in a different form. They are factored into the 
hard part, the collinear part, and the soft kernel which consists of the original soft contri- 
butions and the removed soft part from the collinear part through the zero-bin subtraction. 
Then each term except the PDF is IR finite, and the resummation of large logarithms can 
be performed using the renormalization group equation. We show the factorization theo- 
rems in this picture in deep inelastic scattering (DIS), Drell-Yan (DY) processes near the 
endpoint, and present the one-loop results explicitly. This factorization proof also applies 
in computing resummation near the partonic threshold in the Higgs production, and it 
forms a robust basis in performing a systematic analysis of the threshold resummation. 

The structure of the paper is as follows: In Sec. 2, the features of SCET are described 
briefly and the factorization theorems in DY, DIS processes near the endpoint are outlined. 
In Sec. 2.1, the detailed factorization theorem is presented for DY processes, and the 
factorization theorem for DIS near the endpoint is shown in Sec. 2.2. The factorization 
theorem for the Higgs production near the partonic threshold is presented in Sec. 2.3. In 
Sec. 3, the soft functions for all these processes are computed at one loop in terms of the 
matrix elements of the soft Wilson lines. The soft functions at higher orders clearly show 
that they include IR divergence and mixing between UV and IR divergences, hence not 
physical in themselves. In Sec. 4, we introduce the collinear distribution functions and the 
PDF as the same matrix elements of gauge-invariant collinear operators, but defined at 
different scales depending on the size of the soft momentum involved. Then these functions 
are computed with the corresponding zero-bin subtractions. Here we also explain why the 
need for the zero-bin subtraction in computing the PDF has remained unnoticed in full 
QCD and in SCET. In Sec. 5, the jet function, which describes the collinear particles 
in the final state in DIS, is computed to one loop including the zero-bin subtraction. In 
Sec. 6, all the one- loop results are collected to write the kernel W. And the renormalization 
group behavior of the hard function, the kernel and the PDF is discussed. A conclusion 
is presented in Sec. 7. In appendix A, the computation of the soft function at one loop 
with the nonzero offshellness of the external particles as the IR regulators is presented, and 
the collinear quark distribution functions and the quark PDF in the same regularization 
method is presented in appendix B. And it is also shown that the IR divergence and the 
mixed divergence cancel in the kernel. 

2 Factorization theorems in SCET 

In full QCD, the idea of factorization is probed in detail [1] to separate the hard, collinear 
and soft parts. In describing the soft part, the eikonal approximation is employed to show 
that it is decoupled from the collinear part. In SCET, the procedure of the factorization 
is elaborated since the factorization is achieved at the operator level. The collinear and 
soft parts are expressed in terms of the matrix elements of the operators, and they do 
not interact with each other from the outset by introducing the collinear and soft fields 
without any interaction between them. The scattering cross sections near the endpoint are 
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schematically written as 



da 



H DY (Q,n) ® 5 DY ( £ -,/i) 8 ) f q/Nl (E,n) ® f q / N2 (E, //), 

H D js{Q,fj) ® J(g^/T^,^) ® S-dts(E,h) ® f q/N (E,fi), (2.1) 

#Hig g s(Q,^) ® 5m g gs(^,^) ® f g/Nl {E,n) ® f g/N2 (E,fi), 



in DY [13], DIS [14, 15] and Higgs [16] production processes respectively. H{Q,n) is the 
hard function, Q is the large scale, and u is the renormalization scale. S , (£', ^) is the soft 
function, which is defined in terms of soft Wilson lines. The functions f q /jy(E,[j,) and 
f g w(E, fj,) are the collinear quark and gluon distribution functions defined in terms of the 
gauge-invariant operators at the intermediate scale E ~ Q(l — z) <C Q, where z is the 
Bjorken variable. We reserve the terminology "parton distribution functions" for the same 
matrix elements as those for f q /N and fg/Ni but evaluated at a much lower scale than E. 
The PDF will be denoted as 4> q /N and fyg/N from now on. The final-state jet function 
J appears only in DIS, and describes the collinear final-state particles. The symbol ® 
means an appropriate convolution, and the formulae with the explicit convolution will be 
presented below. 

Eq. (2.1) is the conventional factorization theorem. However, as they stand, the soft 
functions Sjjy, Sbis and Sffiggs contain IR and mixed divergence, so do the collinear func- 
tions f q m 1 , fq/N 2 an d fg/N- The fact that each factorized part can be computed separately 
has a merit, but the existence of IR divergences hinders us from using the evolution of each 
part. It can be cured by reorganizing the collinear and soft parts in such a way that the IR 
divergence is absent except in the PDF. The procedure of reorganization is not arbitrary, 
but is based on the physical principle that consistent separation of the collinear and soft 
modes should be maintained at higher loops. 

Factorization theorems involve disparate scales, and it is more transparent and conve- 
nient to employ SCET to see the physics clearly. There are two scales Q and E ~ Q{1 — z), 
and the effective theories can be constructed by integrating out large scales successively, 
and the matching between the theories can be performed systematically. The first effec- 
tive theory SCETj can be constructed from full QCD by integrating out the degrees of 
freedom of order Q. When we consider DY and DIS processes in SCETj, there appears a 
back-to-back current which can be written in SCET at leading order as 

q^ q = C(Q, rixnyhlYnXn + h.c, (2.2) 

where Xn = Wn£, n is the n-collinear fermion field in SCETi and Xn = WnCfi i s the n- 
collinear fermion. The collinear Wilson line W n is introduced to make the current collinear 
gauge invariant and is given by 

9 



W n = ^2 ex P 



n ■ V + iO 



n ■ A n 



(2.3) 



perm 

where n and n are lightcone vectors satisfying n 2 = n 2 = 0, n • n = 2. The operator n ■ V 
extracts the label momentum, and the bracket means the operator acts only inside the 
bracket. The n-collinear Wilson line Wn is obtained by switching n and n in Eq. (2.3). 
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In the Higgs production, a back-to-back current with gluons is involved, and the detailed 
analysis is given in Sec. 2.3. 

In addition to the UV and IR divergences, there may exist another type of divergence, 
called the rapidity divergence, which occurs when an energetic particle emits collinear glu- 
ons with infinite rapidity. When all the Feynman diagrams are added for a given quantity, 
the rapidity divergences cancel. But they appear in individual Feynman diagrams and a 
regularization is necessary to treat this divergence. It is suggested to include rapidity reg- 
ulators to all the Wilson lines [17]. A simpler method to regulate the rapidity divergence 
is to insert rapidity regulators only in the collinear Wilson lines as [18] 

9 



W n =Y, exp ~ _ ,p , r ~ . n n ■ A n . (2.4) 
V n • V + o n + iO 



perm 



The collinear Wilson line W n can be heuristically derived by attaching n-collinear gluons to 
all the collinear or heavy particles not in the n direction, and integrating out the interme- 
diate states. If these particles are on the mass shell, the gauge-invariant collinear Wilson 
line in Eq. (2.3) is obtained [19]. If those particles are slightly offshell, the collinear Wilson 
line takes the form in Eq. (2.4), where 5 n can be related to the combined offshellness of all 
the collinear and heavy particles not in the n direction. In this case, the gauge invariance 
of the collinear Wilson line is broken. But as long as the sum of all the diagrams is gauge- 
invariant, the introduction of the offshellness can be regarded just as an intermediate step. 
On the other hand, the regulator 5 n can be regarded simply as a regulator with no relation 
to the offshellness. The gauge-invariant formulation with the regulator can be probed, but 
we use Eq. (2.4) as it is to regulate the rapidity divergence. We will show the dependence 
on the rapidity regulator and its cancellation, when summed, at one loop below. 

In Eq. (2.2), the collinear field £ n is redefined to be decoupled from the soft interactions 
as [5] 

£n ^ ^n^ni A n > Y n A n Y n , (2-5) 

where the soft Wilson line Y n is given by 

9 



Y n =Y, ex P ■ A * ■ ( 2 - 6 ) 

^— ' L n ■ 7c + iO 



perm 



Here TV 1 = id^ is a derivative operator on the soft field, which extracts soft momentum. In 
Eq. (2.6), the soft Wilson lines are collectively written as Y n , but they should be determined 
according to the kinematic situations of the scattering process. The detailed prescription 
for the soft Wilson lines is discussed in Ref. [20], and it will be used in defining the soft 
function. 

The matching between QCD and SCETt is performed by comparing the matrix ele- 
ments of the same quantities at the scale Q. The Wilson coefficient C(Q 2 ,/u) in Eq. (2.2), 
for example, is the matching coefficient for the back-to-back current. Since the IR diver- 
gence is cancelled in the matching, the Wilson coefficients are free of IR divergence. Now 
we go further down below E, and the degrees of freedom above E are integrated out to 
produce the final effective theory, SCETn- The matching between SCETt and SCETn can 
also be performed systematically. 
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The key ingredient in adopting the successive effective theories is to relate the collinear 
quark and gluon distribution functions f q /jq and / fl /jv i n SCETi to the quark and gluon 
PDF cf) q/N and <j> g/N in SCET n as [22] 



fg/N(x,v) = 




dz 



dz 



z 



z 



K qq {z,n)4> q / N {x/z,n), 



K gg (z,n)(j) g / N (x/z,iJ,). 



(2.7) 



We distinguish the collinear distribution functions f q /N an d fg/N from the standard PDF 



are evaluated at different energy scales. The collinear distribution functions are evaluated 
at an intermediate energy scale E (Aqcd <£<Q), while the PDF are evaluated at the 
scale much below E, but much larger than Aqcd- As will be shown below, this distinction 
is important. We call K qq and K gg the initial-state quark and gluon jet functions respec- 
tively. The initial-state jet function looks like the matching coefficient for the collinear 
matrix elements between SCETi and SCETn, but strictly speaking it is not. It includes IR 
divergences unlike the Wilson coefficient C(Q 2 , /j,). The initial-state jet function is actually 
the difference in the contributions of the soft modes in the two effective theories. 

The soft function in SCETi describes the contributions of the soft momentum of order 
E ~ Q(l — z), hence the loop momentum of order E ~ Q(l — z) should be subtracted from 
the collinear part to avoid double counting through the zero-bin subtraction. In SCETn, the 
corresponding zero-bin subtraction should be performed with the momentum much smaller 
than Q(l — z). This small momentum is sometimes referred to as the ultrasoft (usoft) 
momentum. These two kinds of zero-bin subtractions are different, and their difference 
yields the initial-state jet function. If there is no distinction between SCETi and SCETn, 
for example, away from the endpoint, the collinear distribution function and the PDF are 
identical and K qq (z,fi) = K gg (z,fi) = 5(1 — z) to all orders in a s , and there are no soft 
functions. 

The initial-state jet function contains IR divergence, so does the soft function. How- 
ever, the sum of these two functions is free of IR divergence. Furthermore the mixing 
of the UV and IR divergences is removed only in this sum. Therefore only the sum is 
physically meaningful and constitutes a component in the factorization theorem. What we 
propose as the proper factorization theorem is to use Eq. (2.7) to express the conventional 
factorization theorem Eq. (2.1) as 



da DY = H DY {Q,V) ® S DY (E,fi) (g> K qq (E,fi) ® K qq (E,/j,) <g> <j> q / Nl (ii) <g> <p q / N2 (fi) 

= #dy(Q, AO ® Wr> Y (E,(x) <g> (f> q/Nl (/i) ® ? /7v 2 (/i), 
dams = H ms (Q,n) <g) J(QVl - z, fj,) (8) Sbis {E, n) <S> K qq (E,fi) (g> 4> q / N {v) 
= H ms {Q^) ® J{Q\Jl - z,fi) 8) Wm S (E,n) <g> <j> q / N (fi), 
do-ffiggs = H mggs (Q,fi) <g> S mg gs(E,n) (g> K gg (E,fi) K gg (E,fi) ® 4> g/Nl (n) <t> g /N 2 (n) 

= H mggs (Q, (j,) <g> Wm gg s{E, m) <8> fl /jVi O) ® 4> g /N 2 (2.8) 



(f) q /N and 4> g /N though the operator definitions for both are the same, but because they 
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for DY, DIS and the Higgs production processes respectively. Here we define the soft 
kernels W as 

W DY (E, fi) = S DY (E, n) ® K qq (E, fj) ® Kgq(E, fj,), 
W ms {E,ri = S Dls (E,ii)®K qq (E,n), 
W mggs (E, fi) = S mggs (E, /i) <g> K gg (E, n) ® K gg (E, (2.9) 

Though the initial-state jet function is included in the kernel W, we call this the soft kernel 
since the initial-state jet function consists of the soft limits of the collinear parts, as will 
be explained below. Eq. (2.8) is our new factorization theorem, in which each factorized 
part can be systematically computed in perturbation theory. Furthermore, all the factors 
except the PDF are IR finite [22]. The explicit convoluted form and the one-loop results 
are shown in this paper. 

One may wonder why we go through this labyrinthine procedure, while we know that 
the PDF in full QCD can be obtained by a straightforward calculation without the problem 
of IR divergence. It will be shown explicitly how full QCD works later. But to put it simply, 
the answer lies in the fact that full QCD computations correspond to the calculation of 
the PDF in SCETn. In SCETn, the usoft zero-bin contribution simply vanishes, hence no 
effect results from the usoft zero-bin subtraction. The details will be explained in Sec 4. 

2.1 Drell-Yan process near the endpoint 

Let us consider the inclusive Drell-Yan process p(Pi)p(p2) 1 + 1~(q) + X(px), where l + l~ 
are a lepton pair and X denotes hadrons in the final state. We define the structure function 
F dy (t) as 

F DY (r) = -N c j ^_e{q Q )5{q 2 - sr) J d A ze~^ (i^l^nO)^^), (2.10) 

with the number of colors N c , and r = Q 2 /s, where Q 2 is the invariant mass squared of the 
lepton pair, and s is the hadronic center-of-mass energy squared. Near the endpoint r — > 1, 
the final-state particles are either soft with the interaction, or n- and n-collinear without 
the interaction. The differential scattering cross section is given as da/dr = (ToFdy(t), 
where oo = 4:7ra 2 Q 2 /(3N c Q 2 ) is the Born cross section for the quark flavor / with the 
electric charge Qf. 

Now we express i ? DY( r ) in SCET by using the current in Eq. (2.2), and the final state 
\X) is decomposed as \X) = \X n )\Xn)\X s ) , the n-, n-collinear and the soft states. The 
momentum g M of the lepton pair is given by q = P\ + P2 —px, where Pi 2 is the momentum 
of the hadron Ni t 2 in the n and n directions respectively. These momenta are given by 
Pf* = n ■ Pin^/2, and Pl^ = n ■ P2n^/2 where s = n ■ P\n ■ Pj. The momentum of the 
final-state particles px can also be decomposed as px = Px„ + Px n + Px s , arid q can be 
written as 

q = Pi +P2 ~ (PX„ +PX n +PX S ) = {Pl-PX n ) + {P2-pXn)-PX s =Pl+P2~PX s , (2.11) 
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where p\^ 2 are the momenta of the incoming partons inside the hadrons N\ i respectively. 
From now on, we express Eq. (2.10) in terms of the partonic variables. First the argument 
in the delta function in the partonic center-of-mass frame can be written as 

q 2 — st = q 2 — Q 2 = (px + p 2 ) 2 - 2p Xs ■ (pi + p 2 ) - Q 2 

= § -2 V sV 2 -Q 2 ~s(l-z-^), (2.12) 

where rj = v ■ px s = (n ■ px s + n ■ px s )/2 = Px s i an< ^ z = Q 2 /^ with the center-of-mass 
energy squared s for the partons. We neglected p x in the first line. Near the endpoint, 
t < z < 1, t — 7-1, and higher powers of 1 — z are neglected. 
The structure function can be written in SCET as 

= -N c H DY (Q,n) J d y^(N 1 N 2 \x n Y^Y nX n 

XXnSfa + - * + 2 ^)^ y ^ 5 (yi - ffj^XnlAW, (2.13) 

where Hdy{Q^) = \Cdy{Q, ^)\ 2 is the hard function and s = yi]j2S. From the last ex- 
pression in Eq. (2.13), the soft interactions are decoupled, and the n- and n-collinear parts 
are also decoupled since they no longer communicate to each other in SCET. The collinear 
matrix elements, after the operators are Fierz transformed and decoupled, can be written 

as 

where a, (3 are color indices, and a, b are Dirac indices. Eq. (2.14) defines the collinear 
quark and antiquark distribution functions for the incoming partons as 

fq/Nifri) = (Ni\x n -5{xin ■ Pi -n-V)xn\Ni), 

f q /N 2 (x 2 ) = (N 2 \x n p(x 2 n -P 2 + n- V^ X n\N 2 ). (2.15) 

These collinear distribution functions are defined in SCETj, that is, the matrix elements 
are evaluated at the scale above E = Q(l — z). 

Finally the factorized form for Fdy{j) is written as 

Fdy{t) = H DY (Q,n) I ——f q /NAyi)fq/N 2 (y2) s VY(z,fJ>) 

J yi V2 y/ q ' 

= H dy (Q,ij)J ^ D y(^)F w -(^), (2.16) 
where the soft function Sdy(z,h) is defined as 

SvY(z,fi) = ^(0\irY^YnS(l - z + Y_ty n |0>, (2.17) 
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and is normalized to 5(1 — z) at tree level. The soft Wilson lines and Y n are chosen 
such that the antiquark and the quark come from — oo [20], and their hermitian conjugates 
are employed in the left-hand side of the delta function. The function F qq is given as 

fJ-) = f —U/NMh/NA—)- (2-18) 
KzJ Jr/z Vl \zyiJ 

Eq. (2.16) is the result of the conventional factorization. 

In our factorization scheme, we further express the collinear function in terms of the 
PDF as 

f 1 dz 

fq/N(x,fj) = / —K qq (z,fl)(j) q/N (x/z,fl), 
J X z 

f 1 dz 

fq/N(x,fj) = / —K qq (z,n)4> q/N (x/z,[i). (2.19) 

J X z 

By shuffling the order of integrations, and changing the integration variables, Fdy is written 

as 

" x dw f 1 dx , , „ , / r 



f dw f dx { 
Fdy(t) = H dy (Q,h) / — / — (f> q j Nl [x^)(f>qj N A 

Jt w Jt/w x ^ 

>j " jK qq (t,v)K qq (^») 



— Sby(z,h) 



Hdy(Q,v) [ —W DY {w,^) [ —4> q/Nl (x,fi)4> q/N ( — ,fi 

Jr W J t/w X H < q < Z \XW 

H DY (Q^) / —W m (w,li)*J L ,lA, (2-20) 



where Wdy and & qq are given as 



- 1 dz 

Wvy(w,(jl) = / — Sby(z,IJl) 



J w/z jK qq (t,v)K qq (^,v), 

= l /w ^/NM^m^^)- ( 2 - 21 ) 

Eq. (2.20) is our new factorization theorem for DY process near the endpoint. The structure 
function is written as the convolution of the hard function, the soft kernel Wdyj an d the 
product of the PDFs. And it will be shown that Wdy is HI finite at one loop. 

2.2 Deep inelastic scattering near the endpoint 

Consider the inclusive deep inelastic scattering e(k\) + p(P) — > e(&2) + X(px) near the 
endpoint. The hadronic tensor W^ u is written as 

W ^ = / d 4 ze^ z (N(P)\J^(z)\X(p x )){X(p x )\r(0)\N(P)) (2.22) 

= ^ E^W* + P ~ Px)(N\jrt\X)(X\r\N), 
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where q = hi — k\ is the momentum transfer from the leptonic system. The summation 
J2x includes all the phase spaces of the final-state particles. Here we consider the elec- 
tromagnetic current only, but the inclusion of the weak current is straightforward. In the 
Breit frame, the incoming proton is in the n direction, and the outgoing particles are in 
the n direction. The analysis near the endpoint in DIS using SCET was first performed 
both in the Breit and the target-rest frames in Ref. [21]. 

The momentum q and the momentum P of the incoming proton can be written as 

77 A* 77 A* 77 A* 

<f = y<3 - y<9, P» = n-P-, (2.23) 



and the Bjorken variable is defined as 

Q 2 Q 



(2.24) 



2P ■ q n-P 

The endpoint corresponds to the limit x ~ 1. 

The hadronic tensor can be expressed in terms of the structure functions as [14] 

= -g^Fi + v"v u F L , (2.25) 
with = (n^ + n^)/2. The structure function F\ is given by 

Fi = -\^W m = -^Y,(2^ 6(4) ( ( l + p -Px)W j \X)(X\J,\N) (2.26) 

x 

= ~r E( 27r ) V4) (3 + p - px)\cmq, ^)i 2 

X 

x j\y(N\x n Yh^nXn\X)(X\x n Yh^ 
where the current J M in SCET is given by 

J" = Cdis(Q, »)XnYh±YnXn- (2-27) 

Here the soft Wilson lines are chosen such that the incoming quark comes from — oo, and 
the outgoing quark goes to oo. 

We now define the two collinear functions in the n and n directions, which do not 
communicate to each other. The final-state jet function in the n direction is defined as 

^XnlXnHXnlxn = \J ^Jfr*.). (2.28) 

Xn 

The jet function is a function of n ■ px n only. The collinear distribution function in the n 
direction is defined, as in Eq. (2.14), by 

Y / (N\(xX\Xn)(X n \6(y- 1 ^)(xnf b \N) = ^f^^JiMv)- (2-29) 
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In terms of these two collinear functions, the structure function F\ can be written as 

Fx(x) = \C ms (Q^)\ 2 C-f q /N{y^) 



a y 

rQ(i-x/y) dn . k / x n-k \ 

x / — JS'i,|i )S ms 1 pr-./x , 2.30 

Jo 2tt V y Q J 

where the soft function Sdis is defined as 

5dis(1 ~z) = ^(0|tr[Y n ty fi <y(l - z + ^)y_ty n j | 0) . (2 . 31 ) 
Changing the variables n ■ k = (1 — w)Q, and defining 

J(w,/jl) = j-j(Q(l -w),fij, S D is(—,ixj=wSms(w-z), (2.32) 

The jet function J(w,fi) and the soft function Sj^isiw, fi) are normalized to 5(1 — w) 
respectively at tree level. Eq. (2.30) is written as 

f ^ dz / x \ dw / z \ 

Fi(x) = H ms (Q,Li) — f /N / — J(w,ij,)Sbis — ,M , ( 2 -33) 

J x z ' \z / J z w \w / 

where i?Dis(Q) A 4 ) = |Cdis(Q> A 4 )! 2 ; an d subleading terms of Sbis m w is neglected near the 
endpoint. 

Eq. (2.33) is the result of the conventional factorization theorem near the endpoint. In 
our formulation, as was done in DY process, we go further by invoking Eq. (2.19) to write 
Eq. (2.33) as 

*i(x) = flbis(9,M) f - t t —JHSbis(-). (2.34) 

Jx z J x/z y \yzJ J z w \wJ 

Changing the order of integrations successively, and redefining the integration variables, 
the factorized structure function is given as 

f dw f dv r f du / v \ i / x 

F 1 (x) = H(Q,n) — J(w) — / —S D is{u,iJ,)K qq [-,fj,)(j) q / N [ — ,/i 
J x w J x/w v ij v u HH \u J\ Hi \vw 

f dw f dv / x \ 

= H(Q,n) —J{w) —Wvn(v,ri<l> q/N (—iA, (2.35) 

Jx W J x/w V \VW J 

where WdisC^jAO is given by 

W ms (v,fj,)= —S Dls (u,v)K qq (-,v). (2.36) 

Eq. (2.35) is our new result of the factorization theorem in deep inelastic scattering. The 
soft kernel Wdis is IR finite, and surprisingly enough its radiative corrections to all orders 
vanish and Wdis( v ) = ^(1 — v). The argument about why it is true to all orders in a s is 
given after the one-loop corrections are presented. 
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2.3 Higgs production near the partonic threshold 

The factorization proof for the Higgs production is similar to that for DY process except 
that the initial-state particles are gluons instead of a qq pair. Therefore this process involves 
the gluon PDF and the initial-state gluon jet function. The effective Lagrangian for the 
Higgs production after integrating out the top quark loop is given by 

C eS = C t (m t ,fi)^-G a ^G^ a , (2.37) 

where v is the electroweak vacuum expectation value, mt is the top quark mass, H is the 
Higgs field and G°^ v is the field strength tensor for gluons. The coefficient Ct(mf,/x) to 
second order in a s , and in the heavy top quark mass limit, is given by [23, 24] 

C t (m t , M ) = w (l + — -y (2.38) 

The gluon field strength tensor is matched onto SCET as 

G%G^ a -> -2C H (Q,»)B^yty n B£ x = -2C H (Q,»)O g , (2.39) 

after decoupling the soft interaction. In this process, Q = mji, the Higgs mass. Here £>n M 
is defined as 

B^ = i[n • VWliD^Wn) = in a g%WlGfW n = in a g^T a (Wl) ab Gf' b . (2.40) 

In the final expression, the collinear Wilson line W n is written in the adjoint representation 
rather than in the fundamental representation. That means that the generator T a is given 
by (T a ) bc = -if abc . The collinear Wilson lines can be expressed either in the fundamental 
representation or in the adjoint representation. Both approaches are equivalent, but the 
adjoint representation is employed here since it makes the expression in the factorization 
formula look similar to that in DY processes. The Wilson coefficient Ch is the matching 
coefficient between the full theory and SCET and is given to order a s by [25] 

Ch(Q, A*) = 1 + (- In 2 i£ + \^ - Hit In ^) , (2.41) 

where Ca = N c . Then the effective Lagrangian in SCET is given by 

£scet = -^Ct(m t ,^C H (Q,^HO g = -C{n)HO g . (2.42) 

The kinematics of the Higgs production process p(P\)p(P 2 ) — > HX is similar to that 
in DY process. The momentum q of the Higgs particle is given by 

q = Pi + P2-PX =Pi+P2~Px s , (2.43) 
as in Eq. (2.11). The cross section for the Higgs production process is written as 

a(pp -+ HX) = - I d\5^(P 1 + P 2 -q- p X )S(q 2 - Q 2 )%°) 
s x 

x\C{^)\ 2 {N x N 2 \0\\X){X\O g \N x N 2 ) 
^-ic^i 2 I dx 1 dx 2 (N 1 N 2 \B^ a yl ab y^ls(l - 2,1 



X 



±ud 



|7V 1 iV 2 ),(2.44) 



- 13 - 



where rj = v ■ px = Px an d z = Q 2 /s (r = Q 2 /s ) with the center-of-mass energy squared 
for the partons (the hadrons). The partonic threshold corresponds to the limit z — > 1. As 
in DY process, the soft interactions are decoupled, so are the n- and n-collinear parts. The 
collinear matrix elements are written as 



S(xi 



n-V 



|iVl) = 2(^-l) x( "- Pl)2/ ^ (x) - 



Inverting Eq. (2.45), the collinear gluon distribution amplitude / fl /jVi is written as 

n-V 



fg/NA x ) 



n-Pi 



\Ni). 



(2.45) 



(2.46) 



For f g /N 2 i n an d n are switched in Eq. (2.46). 

In terms of the collinear gluon distribution functions, the cross section is written as 

<r(pp^HX) = a J dx 1 dx 2 H mggs (Q, ^)S mggs (z, n)f g / Nl (xi, n)fg/N 2 (x 2 , (2.47) 

where the soft function Sm gg s(z, n) is defined as 

Sm ggs (z,u) = ^^<0|trXto(l " z + 2 ^)yty n \0). (2.48) 

The Born cross section <to is given by 

2Tr\C t (m t ^)\ 2 



O"0 



v 2 {N?-l) ' 



(2.49) 



and H mggs (Q,u) = \C H (Q,u)\ 2 . 

Here we invoke the relation Eq. (2.7) to write the cross section as 

a(pp -> HX)/a = H mggs (Q,fi) j —f dx<p g/Nl (x, fJ,)<f> g /N 2 (— , (jl 

Jt ^ Jt/w \XW 

^S mggS (z,fi) [ dtKgg(t,fl)Kgg(^- H 



Hm ggs (Q,fJ-) / — W mgg s(w,n) / dx(/) g / Nl (x, fj,)(f) g/N J — 

Jt ^ Jt/w \XW / 



W 



(2.50) 



where Wniggs and $»» are given as 



W mgg s(w,n) = ( — Smggs(z,n) ( dtK g g(t,(i)K gg (^,[i\, 

Jw z Jw/z ^ z t ' 



w 



i 

t/w 



dX<j)g/ Nl (X, fl)4>g/N 2 [ — •/' 



\xw 



(2.51) 



Eq. (2.50) is our new factorization theorem for the Higgs production near the partonic 
threshold. 
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3 Soft function 



In the conventional factorization scheme, the most problematic part is the soft function 
because it includes the IR and mixed divergences. Unless these divergences are removed, 
the soft function is not physical and the evolution via the renormalization group equation 
is meaningless. Especially the real gluon emission in DY process has only IR divergences. 
From this section, the ingredients in the factorization formulae Eqs. (2.20), and (2.35) 
are computed at one loop. The UV and IR divergences appear unequivocally in the soft 
functions, and we compute the one- loop corrections to the soft functions. 
In DY process, the soft function is given by 

Svy(z,h) = ^-(0|trY n tY^(l - z + 2 ^y^Y n \0), (3.1) 

The Feynman diagrams for the radiative corrections of the soft function at one loop are 
given in Fig. 1. Fig. 1 (a) with its mirror image corresponds to the virtual correction, 
and Fig. 1(b) is the real gluon emission. The computation is performed by employing the 
dimensional regularization to regulate both the UV and IR divergences with the spacetime 
dimension D = 4 — 2e and the MS scheme. The results with another regularization scheme, 
in which the dimensional regularization is employed for the UV divergence, and the IR 
divergence is regulated by the offshellness, are presented in Appendix A. 
The corresponding matrix elements are given as 

n - 2^ /^ 2 e 7E v f d D l 1 a s C F / 1 1 \2 

i 2 e"' E \£ f d D l 1 „,.,, / 2vl 



ol s Cf 



2tt 



*(!-*)(-*- + — ln^+21n 2 ^-^ 
v em e iR Q Q 4 



2 /I , M\ ,/ln(l 
— — + 21n^ +4(^- 



l-z) + \ei R QJ \ 1 — z J + 



(3.2) 



In this computation, the function 1/(1 — z) 1+t appears. It diverges at z = 1, and the 
divergence is definitely of the IR origin. In terms of the plus distribution functions, it can 




Figure 1. Feynman diagrams for soft functions at one loop (a) virtual corrections and (b) real 
gluon emission. 



be expanded in powers of e as 



1 l r/ 1 /ln(l-z)\ , . 

-£(1-*)+ . , -e +■■-. (3.3) 



(l-z)^ e v ' (1 

Here M* DY is the result with y n (Y^) in the left- (right-) hand side in Fig. 1 (b), and the 
result with and Yn should be included. They are hermitian conjugates to each other. 
The virtual correction has the UV and IR divergences, and the mixed divergence. Note 
that there are only IR divergences in the real gluon emission. 

The total soft contribution for DY process at one loop is given as 



4 1 Y (z) = 2Re(M s V + M s fe DY ) 



a s C F 



7T 



, 2 1 , fl 2 1 , 2 /i 2 7T 2 

6(1 -z)( — 2— + + — In % + - In 2 % - - 

K £uv £ uvem e m Q 2 2 Q 2 4 

2, ,wi z] 



(l-zj+VeiR Q 2 / V l-z 



(3.4) 



In Eq. (3.4), the finite terms are the same as the result in Ref. [26]. But the divergent 
terms are mixture of IR and UV divergences contrary to the argument in Ref. [26] , in which 
they claim that the IR divergences are cancelled due to the KLN theorem. Our explicit 
calculation shows that it is not true. The KLN theorem holds when there is no kinematic 
constraint on the real gluon emission. However, as can be seen in M^ DY for the real gluon 
emission, there is a constraint for the soft momentum of the real gluon specified by the 
delta function, which is different from the virtual corrections. Therefore the cancellation 
of the IR divergence is bound to be incomplete. 

The soft function for the Higgs production Smggs^) is defined as 

Sm sgs (z) = ]N ^ I (0|tr^5(l - z + 2 ^)^n|0), (3.5) 

where Q = mn, and the soft Wilson lines y n , and 3^n are in the adjoint representation. 
•Sffiggs^) is normalized to 8(1 — z) at tree level. In this expression, Eq. (3.5) resembles 
Eq. (3.1), and the only difference is the representation of the soft Wilson lines. This means 
that the radiative corrections for <Sffiggs is the same as Sdy except the color factors. By 
explicit computation, it is true and the corresponding radiative corrections for 5*Higgs ar e 
obtained from Eq. (3.2) by replacing Cf by Ca- Accordingly, the one- loop correction for 
Sh iggs is given by 

(i) ,._a s C A \, ( , ,f 1 , 2 1 a 2 1, 2 M 2 ^ 



SffiggsW = — f(l-*)f--a-+ + — ln^ + -In' 

vr L V efj V euvem em Q l 2 Q 2 



2 / 1 i [J* 2 \ „/ln(l 
— — + m^r +4 4- 



(l-z)+Ve m Q 2 J V 1 
The soft function Sbis m DIS near the endpoint is given by 



(3.6) 



Sdis(*) = ^(0|tr[y n ty ft< yfl - z + r ^)Y^Y n ] |0). (3.7) 



Nr. 



Q 
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The Feynman diagrams for the one- loop correction is shown in Fig. 1, except that Yn is 
replaced by Y^. The matrix elements are given as 



M, 



s.DIS 



2vr Veuv em/ 

at s C F / 1 
2vr V euv em ^ 



1 MY 1 1 i » 



(1-*)- 



(3.8) 



Unlike DY process, the real gluon emission contains UV divergence. The total soft contri- 
bution in DIS at one loop is given by 



Sg} s (z) = 2Re(M» DIS + M s & >DIS ) 



(3.9) 



—(---)[(--- 2 ln i) ^ - *> + jy 1 - 1 

vr Veuv eiR/ LV ^UV 2 (1 - z). 



All the soft functions contain IR, UV divergences as well as mixed divergence. There- 
fore they are not physically meaningful as they are. As claimed, only after the initial-state 
jet function is added, the IR and mixed divergences disappear. 

4 Collinear distribution functions and PDF 

4.1 Collinear quark distribution function and PDF 

The radiative corrections of the collinear quark distribution function, defined as 

f q/N (x, fi) = (N\x n ^S(n -Px-n- P) Xn \N), (4.1) 

will be computed explicitly at one loop. The Feynman diagrams for the collinear function 
at one loop are shown in Fig. 2. Fig. 2 (a) is the virtual correction, and Fig. 2 (b) and (c) 
are real gluon emissions. The mirror images of Fig. 2 (a) and (b) are omitted, which are 
given by the hermitian conjugates of the original diagrams. 




Figure 2. Feynman diagrams for collinear functions and PDF at one loop (a) virtual corrections, 
(b) and (c) real gluon emission. The mirror images of (a) and (b) are omitted. 



The matrix elements are written as 

'l^e^y.. f d D l n-(p-l) 



M a = 2ig 2 C F (^y5(l-x) J 



{2tt) d P(l-p) 2 (n-l + h) 



a 2^ //i 2 e 7E \ e f d D l n-(p-l) / n-l\ c .^ 

"-VC^V-^/^^^.-gift (4.2) 

These are computed using the dimensional regularization for both UV and IR divergences 
with p 2 = in the MS scheme. The results with the dimensional regularization for the UV 
divergence, and the nonzero p 2 for the IR divergence are presented in Appendix B. 
The results of the computation are given by 

^*(l-z)(J--J-Vl + ln * 



2ir " Veuv eiR/V n-p 



M a 

5(1 _ x)ln= ^ + 

2tt \e uv eiR^ V n-p (1 

*i = ?^(l-,)(J— ±). (4.3) 
27r v euv eiR y 

The rapidity regulator 5\ is employed in the computation. As can be seen, M a and Mf, 
depend on this regulator 5\, but this dependence is cancelled in the sum M a + Mf, even 
without the zero-bin subtraction. This is to be contrasted with the transverse-momentum- 
dependent collinear function, where the cancellation of b\ is achieved only after the zero-bin 
subtraction [18]. 

In computing the zero-bin contributions, we neglect all the components of the loop 
momentum I compared to n ■ p. From Eq. (4.2), they are given as 

i 2 eT E \e , s f d D l 1 



Mr = -2^C F (^)V-)/ 



(2vr) D l 2 n ■ l(n ■ I + 5i) 



2vr Veuv £ir/ Veuv A* 

= wcJif-X f f' „ \ ^ 5(l - x - ^ W 

b ' s V 4-7T / 7 (2n) D n ■ l(n ■ I + 5i) V n-p) 



^(—--)(-S(l-x)lnA + 
2ir Veuv £tr / V n-p 



1 



Wv eiR/V n-p (1-x). 

= 0. (4.4) 

The zero-bin contribution is subleading and is neglected. Note that the distinction 

between the soft and usoft zero-bin contributions appears in M 6 (0 S } . In SCETi, the soft 
momentum is of order Q(l — x), which is of the same order as the loop momentum n ■ I 
in the zero-bin contribution. However, there is no distinction between the soft and usoft 
contributions in because the integral remains the same irrespective of the size of the 
loop momentum. can also be neglected in the usoft limit. Therefore the soft and 

usoft zero-bin contributions are different only in Mj,, which corresponds to the real gluon 
emission with soft or usoft momentum. 
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M a 


= M a 


- Mf = 


M b 


= M b 


- M<°> = 

b,s 


M c 


= M c 


Q! s Cf , 

2tt 1 



The collinear part with the soft zero-bin subtraction is written as 

^ (1 _ x) (J__J_)(J_ + 1 + ln JL 

2vr Veuv £ffi7 Wv n-p 

a s C F t 1 1_\ 

2tt V euv em ' ' 



l-x) . (4.5) 

v euv em ' 

Note that various combinations are independent of the rapidity regulator Si though indi- 
vidual diagrams depend on it. As was first noted, the sum of naive collinear contributions 
is independent of <5i. So is the sum of soft zero-bin contributions, hence the true collinear 
contribution with the zero-bin subtraction. It is also true that the usoft zero-bin contri- 
bution is also independent of Si, in fact, it vanishes. It turns out that the same result is 
obtained without introducing Si in the beginning. But Si is included in the calculation 
to show how the cancellation occurs explicitly, and this method has been used also in 
calculating the transverse-momentum-dependent collinear distribution function [18]. 
The collinear distribution function at one loop is given by 

/%(*, = 2Re {M a + M b ) + M C - ^ (- -) S{1 - x) 

qi 47r Veuv eiR y 

_2^(J__J-)r(J- + « +2h jL) <(1 _ l) _ 1 _ I i (4.6, 

27r Veuv em/LVeuv 2 n-p) J 

where the last term in the first line comes from the self-energy of the fermion field £. 
The radiative correction contains the UV and IR divergences, furthermore there is 
the mixing of the UV and IR divergences. Due to the IR and mixed divergences, it is not 
appropriate to discuss the evolution of the collinear distribution function. Our factorization 
procedure removes the mixing of the UV and IR divergences as will be shown later. 

In SCETn, the zero-bin subtraction must be usoft, and the only difference from the 
soft zero-bin subtraction appears in M b . The usoft zero-bin contribution for M. ^ is given 
by 

w(0) , 2^ ^ 2 e 7E \% /1 v f d D l S{1 2 ) 

= °^8{l- X )(±- --)(— -In^' 
2vr Veuv e m ) Veuv M 

and the corresponding collinear contribution is given by 
r ( ) = a s C F / 1 1_\ 



(4.7) 



M b ,» s = M b -M b 

2ir V euv eiR ' 



8(l-x)[-—-]nJ-) + 7 -^- T - . (4.8) 
euv n-p) {l-x). 



Note that M aiUS = M ajS = — M bjXIS , hence M QiUS + M b)US = 0, meaning that there is no usoft 
zero-bin contributions in the PDF. 

The PDF <j& N (x,fi) at one loop is given by 

<jP-hx, m) = 2Re {M a + Mb us) + M c - (— -)s{l - x) 

q/ 47r Veuv effiy 

+ (4.!-) 
27r Veuv e m )\2 (l-x)+J 
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Unlike f^ N , the radiative correction (f> q l ) N does not involve the mixed divergence. The 
result is the same as in full QCD. The IR divergence is absorbed in the nonperturbative 
part of 4> q /N, and the UV divergent part yields the anomalous dimension of the PDF, which 
governs the renormalization group behavior. 

The fact that the usoft zero-bin contributions in the PDF vanish is responsible for why 
full QCD results can be obtained from naive collinear contribution only. Of course, the 
full QCD results hold away from the endpoint region, and near the endpoint. It can be 
clearly explained in SCET. Away from the endpoint region, there is no intermediate scale 
and there is no need to construct effective theories successively. And there are no delta 
functions in the soft functions in Eqs. (3.1) and (3.7). Then the soft functions are just 
one, meaning that there are no soft contributions to all orders in a s . And the radiative 
corrections of the collinear part at low energy scale are exactly those of the PDF given by 
Eq. (4.9). However, there is no way to avoid the steps described above near the endpoint 
region. 

The relation between the collinear function / 9 /jv an d the PDF g /jy is given by 



/ ' dz / x \ 

fq/N(x,fJ,)= —K qq (z,fl)(j) q/N l-,fl), (4.10) 
J rp Z \ Z / 



where K qq (z,n) is the difference between the zero-bin subtractions in the collinear matrix 
element between SCETi and SCETn. At one loop, it is given by 

JfgW) = 2Re(-M<J + M<°>) (4.11) 



L CJ(J Vlt*) 6,8 ' 6,us 

^(J--J^fl-^J-+ln^ 
7T V euv eiR 



5(1 -z){ + ln_ 

V euv n ■ p 



We stress again, K qq is not the matching coefficient between SCETi and SCETn- It 
contains IR divergence, and the mixed divergence too. There also appear IR and mixed 
divergences in the soft function. Combining these two, there will be neither IR nor mixed 
divergence, as is explicitly shown at one loop. 



4.2 Collinear gluon distribution function and PDF 

The collinear gluon distribution function is defined as 

f g /N(x,ii) = -^—(N(P)\B^ a 5(xn- P -n-V)B%\N(p)), (4.12) 

which is normalized to f g /Ni x ) = S{1 — x) st tree level. The Feynman diagrams for the 
collinear gluon distribution function at one loop is shown in Fig. 3. 
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Figure 3. Feynman diagrams for collinear gluon distribution functions and PDF at one loop (a) 
virtual corrections, (b) and (c) real gluon emission. The mirror images of (a) and (b) are omitted. 



The matrix elements in the background gauge are given as 



M a = 2ig 2 C A n ■ P 5(l - x) 



4ir 



d D l 



(2tt) d l 2 {l + p) 2 Vn-l -<5i n-(l + p) + S 1 



.2/ 



In 



Si 



*(l-s)(J-J_) 
v euv eiR/ 

[n • (I — p)] 2 rn • (p — 21) 2n • p 



(l-pf 
a s C A { 1 



n ■ p 



n-l + Si 



M r . = 



2tt V euv em 



1 — x) In 



Si 



+ 



• , + x(l - x) - ^ ) , 

re • p (1 — x)+ 2/ 

^(^)7I>M»-'-<'-*H<^(t»-'"-'- i 3 



Q s Ca / 1 M 1 - 

Wv e IR /\2 x 



(4.13) 



The soft zero-bin contributions are given as 
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The usoft zero-bin contribution only differs in the real gluon emission of Mb, and it is given 
as 



(0) = OsCa _ x) / J_ _ M / J_ _ M 
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Therefore the matrix elements of the collinear gluon distribution function at one loop with 
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the soft zero-bin subtraction are given by 
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and the usoft zero-bin subtraction M\, )US is given by 
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Combining all the ingredients, the collinear gluon distribution function at one loop is 
written as 
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where the last term in the first line is the self-energy correction of the gluon field in the 
background gauge, and /3q is given by 
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The gluon PDF (fig/N is obtained by replacing Mb,s with Mb jUS in Eq. (4.18), and is given 
by 
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The one-loop correction to <P^ N is the same as the result in full QCD for the same reason 
as in the case of the quark PDF. The initial-state jet function at one loop is given by 



= 2Re(-M 6 ° s + M° J 
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(4.21) 
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Compared to the initial-jet function with quarks in Eq. (4.11), K gg is the same except the 
color factor, satisfying K^q(x)/CF = K g l g (x)/CA- 
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5 Final-state jet function 



The final-state jet function for DIS in the n direction is defined as 

\- \Y \/TT \- ^ f d4pX ' 

2^Xn\Xn){X n \xn = 2 / 



(27T)' 
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(5.1) 



The jet function J(n-px) is a function oin-px only. As defined in Eq. (2.32), the radiative 
corrections to J(z,fi) = QJ(Q(1 — z))/(2n) will be computed. It has been computed to 
two- loop order [27], but we will explicitly present the calculation at one loop to show how 
it can be computed like the collinear distribution function. However, the final states are 
not on the mass shell p\ ~ Q 2 (l — z), the IR divergence is regulated by the offshellness 
p 2 x . In the intermediate calculation, the IR poles appear, but in the final sum they cancel. 

The Feynman diagrams for the final-state jet function at one loop are shown in Fig. 4. 
Fig. 4 (a) and its mirror image are the virtual corrections, and Fig. 4 (b) and (c) are the 
real gluon emissions. The virtual correction of the fermion self energy is omitted in the 
figure, but is added separately in the final calculation. Their matrix elements are given as 
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n-p J (1 — z)- 

To be consistent with the idea of the zero-bin subtraction, the soft contributions should 
be subtracted from the naive collinear calculation. The zero-bin contributions in this case 
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wl 
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Figure 4. Feynman diagrams for the final-state jet functions in DIS at one loop (a) virtual 
corrections, (b) and (c) real gluon emission. The mirror images of (a) and (b) are omitted. 
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are easy to deduce, and the soft zero-bin subtraction is the appropriate procedure. The 
zero-bin contributions from Eq. (5.2) are given as 
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where Mc here is again subleading and is set to zero. The final result for the final-state 
jet function with the zero-bin subtraction is written as 
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Therefore the final-state jet function in DIS at one loop is written as 
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where the last term in the first line comes from the virtual correction with the self-energy 
for the fermion field. In Eq. (5.5), n ■ p is replaced by Q, appropriate in DIS. This result is 
the same as the one in Ref. [15]. It is confirmed that the final-state jet function has only 
UV divergences after the zero-bin subtraction. 



6 The soft kernels W and the renormalization group equation 

Combining all the one-loop results, the factorization formulae in Eqs. (2.8) and (2.9) can 
be explicitly presented. The kernel W to one loop is obtained by plugging all the one-loop 
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results in Eqs. (2.21), (2.36) and (2.51). They are given as 
W DY (z,ri = 8(1 -z) 
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As claimed, the soft kernels are free of IR divergences. These kernels can be computed 
systematically using perturbation theory, and the anomalous dimensions can be obtained 
from Eq. (6.1) to describe the scaling behavior. 

Remarkably the kernel for DIS becomes 6(1 — z) to one loop since 5^ ig + Kqq = 0. 
Therefore the inclusive DIS cross section near the endpoint consists of the hard part, the 
final-state jet function and the PDF <fi q /N at the scale [i < E. It means that the soft 
contribution and the initial-state jet function cancel. Therefore the PDF satisfies the 
ordinary evolution equation, as in full QCD. If we do not combine the soft function and 
the initial-state jet function as was done in conventional factorization approach, the soft 
function and the collinear part f q m = K qq 4> q /N include IR and mixed divergences. 
Then neither the soft function nor the collinear distribution function is physical, and it is 
meaningless to consider the evolution of the collinear distribution function. The solution 
to this problem is our approach. That is, we combine the soft function and the initial-state 
jet function to make an IR-finite kernel, which has no radiative corrections for DIS, and 
the remaining collinear part is the PDF 4> q /N gives the same result as in full QCD, given 
by Eq. (4.9). The difficulty in treating the evolution of the collinear part without the 
reorganization has been discussed in Refs. [14, 15, 28]. 

The fact that VFdis = 5(1 — z) is true to all orders in a s , and it can be shown by a 
simple argument. The soft zero-bin contribution from f q /N 1S obtained by integrating out 
the momenta of order QX. This amounts to attaching a soft gluon to a collinear fermion, 
and making a loop with the soft gluon. This is exactly the procedure to obtain the eikonal 
form of the soft Wilson line and to calculate its loop correction. That is, the soft zero-bin 
contribution is the same as the soft function in DIS. On the other hand, the usoft zero-bin 
contributions vanish in <f> q /j^ to all orders in a s . It is explicitly verified here at one loop, but 
if we look at Eqs. (3.1) and (3.7), the derivative term in the delta function is much smaller 
than (1 — z), hence can be neglected. Then the delta function can be pulled out, and the 
remaining usoft Wilson lines cancel. The usoft function becomes 5(1 — z) to all orders in 
ol s . Therefore in the absence of the usoft contributions, the initial-state jet function, which 
is the negative value of the soft zero-bin contributions, always cancels the soft function to 
all orders in a s in DIS. 

In contrast, it is different in DY processes since the soft function involves an interaction 
between n and n-collinear fermions, while the collinear part which interacts only within 
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each collinear sector does not produce the same soft interaction in the zero-bin limit. 
Note, however, that Wjyy is still IR finite though the soft function involves the interaction 
between different collinear parts, and the initial-state jet function includes the interaction 
only in each collinear sector. The disparity between the soft function and the initial-state 
jet function becomes acute in multijet processes, and it will be interesting to see if the 
corresponding kernel will still remain IR finite in a more general case with multijets. 

The factorized forms of the structure functions involving the kernels W have been 
already shown in Eqs. (2.20) and (2.35). Each factorized function has a nontrivial UV 
behavior, however when we combine all together, the structure functions and the scattering 
cross section should have no scale dependence. 

^ %(T)=0 ' ^ Fi( " ) = °' V m ^ = °- (6 - 2) 

Since all the elements in the factorization theorem are computed to one loop, the evolution 
of each quantity can be derived to next-to-leading logarithm accuracy. 
The Wilson coefficients Cdis(Q 2 ,m) an d Ch{Q,Ij) are given by 
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Using the relation Cdy(Q 2 )/-0 = Cdis( — Q 2 ,fi), and from Eq. (2.41), the hard coefficients 
to one loop are given by 

7T 2 X 

+ t)< 

7tT 2 



Hnis(Q,v) 


= 1 + 


a s CF 1 
2vr ' 


' , 2 M 2 
" ln Q 2 






= 1 + 


OI s Cf 1 

2vr ' 


-ln 2 ^ 

v q 2 


- 3l1 ^ 


Hm g gs(Q, m) 


= 1 + 


0( s Ca 1 
2tt ' 


' 1 2 /" 2 

" ln Q 2 


+vy 

6 / 



(6.4) 



7//dy (A*) = ^DY = ( 2 ln 7p + 3 ) = Tffws , 
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From Eq. (6.4), the anomalous dimensions of the hard functions in DY, DIS and the 
Higgs production processes at one loop are given by 

d a s C F t /x 2 

^ Hdy = - — ( 21n Q 2 
d 2a s C A , l± 2 

^ Hm ^~ — ^ ln Q 2 

In the Higgs production, there is also \Ct(mt, n)\ 2 in ao, which is proportional to a 2 . To 
one loop, the anomalous dimension for \Ct(mt, /u)| 2 is given by 

din |C t | 2 a s C A /ll 2n/\ 

7C < = = -— It - wj ' (6 - 6) 

which comes from the running of a s at one loop 
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For the soft kernel Wdy> Wffiggsj and the PDF (p q /N, 4>g/Ni they satisfy the renormal- 
ization group equation 



d 



dz 



d , f 1 dz 



— 1g{z,n)<j) a / N (-,H 



The anomalous dimensions 7vy(z,^) and 7g )fl (z,jti) are schematically given as 
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where Z is the counterterm and 7 is the anomalous dimension for W or <f> q /N, (j)g/N- As 
discussed, Wdis does not evolve. 

From Eq. (6.1), the anomalous dimension of the kernel, Wdy> is given to one loop by 
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and the anomalous dimension 7w ffl ggs f° r ^Higgs is obtained from yyy by replacing Cf by 
Ca- From Eq. (4.9), the anomalous dimensions for the quark PDF 4> q /^ and the gluon 
PDF are given by 
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where and P gg are the splitting kernels appearing in the DGLAP evolution equa- 

tions for the PDF. The anomalous dimension for the final-state jet function is computed 
as 
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For each process, the sums of the anomalous dimensions near the endpoint x — > 1 are 
given as 



T/fcis'H 1 ~x) + j q (x) + 7j(x) = 0, 
7# DY (5(1 - x) + 7w DY (x) + 2-f q (x) = 0, 
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Combining all these anomalous dimensions, we can see explicitly that Eq. (6.2) holds true 
to one loop near the endpoint. 
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7 Conclusion 



The conventional factorization theorems have been successful in the sense that the effects of 
strong interactions at various stages have been satisfactorily separated to express scattering 
cross sections as convolutions of the high-energy part, the collinear and the soft parts. 
But there has remained a problem since the divergence structure is so intricate that the 
collinear and the soft parts still contain UV, IR, and mixed divergences. Now we have 
better understanding of the origins of the divergences, and this problem is taken care of 
by separating the soft modes consistently in the collinear part. As a result, the factorized 
parts no longer involve problematic IR or mixed divergences. With our new factorization 
scheme, factorization theorems have gained stronger grounds for the theoretical description 
of high-energy scattering. 

Our factorization formula starts with a physical idea that the soft and collinear modes 
should be separated at higher loops, and employs the zero-bin subtraction to realize this 
idea. Our new factorization theorem emphasizes consistent separation of each mode at 
higher loops. In loop calculations, the soft contribution always encroaches on the collinear 
sector since the collinear loop momentum covers the soft region. To ensure the separation 
of the collinear and soft parts, the soft contribution in the collinear sector should be con- 
sistently removed from the collinear sector. Otherwise the overlap is bound to occur in 
every collinear loop calculation. 

The consistent treatment of the zero-bin subtraction results in the appropriate diver- 
gence behavior. Without including the zero-bin contributions from the collinear part, the 
soft contribution contains not only the IR divergence but also the mixing of the IR and UV 
divergences. The mixed divergence is especially troublesome and it should be absent for 
the soft function to have physical meaning. The inclusion of the initial-state jet functions 
K qq or K gg removes this mixed divergence. In addition, it also changes the IR divergence 
to the UV divergences, and the kernels W are physically meaningful. If we naively put 
e = euv = eiR an d identify the poles in e as the UV poles, the soft function Sj^y m Eq. (3.4) 
is identical to the kernel Wjyy in Eq. (6.1). It would be a good mnemonic to identify the 
UV divergence, but physically it does not make sense. Since we know the origin of the 
divergences, the UV and IR divergences can be systematically identified, and a physical 
quantity should be free of IR divergence. 

Let us finally summarize the recipe for our factorization procedure. First, we write 
down the scattering cross section. Second, as in the conventional approach, it can be 
factorized into the hard, collinear and soft parts. In full QCD, they can be scattering 
amplitudes. In SCET, the hard part is obtained from the Wilson coefficient, the collinear 
and the soft parts are defined as the matrix elements of the relevant operators. Third, 
radiative corrections are computed. If there is an intermediate scale, the collinear distribu- 
tion function with the soft zero-bin subtraction above the scale can be related to the PDF 
with the usoft zero-bin subtraction below the scale. The relation is expressed in terms of 
the initial-state jet function. Finally, we combine the soft function and the initial-state jet 
function to yield the soft kernel, which is IR finite. 

Our factorization scheme gives a consistent field theoretic treatment of the UV and IR 
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divergences. This scheme has also been successfully applied to DY processes with small 
transverse momentum [18]. In this case the size of the transverse momentum offers the 
intermediate scale which distinguishes SCETi and SCETn- The initial-state jet function 
takes the form 

fq/N{x,k±,n) = — Kg q Qs._ L ,n)<p q / N {-,{!), (7.1) 

J x Z \ Z / 

where f q /]y(x, kj_, jj) is the transverse-momentum-dependent collinear distribution func- 
tion, and Kj g (k.±, fi) is the initial-state transverse-momentum-dependent jet function. 
Combining the initial-state jet function with the transverse-momentum-dependent soft 
function also yields an IR-finite soft kernel. Therefore our factorization formalism can 
be applied to various high-energy processes. It remains to be seen whether this can be a 
general formalism for factorization proof. A research in this direction is in progress. 

A Soft functions with offshellness 



In the appendices, the soft functions and the collinear distribution function are computed 
in the regularization scheme with the offshellness of the external particles for the IR di- 
vergence. By computing the soft kernel and the PDF in this scheme, we also show the 
cancellation of the IR and mixed divergences. 
For DIS, the soft function is defined as 
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But here instead of Eq. (2.6) for the soft Wilson lines, the offshellness of the collinear 
particle from which the soft Wilson lines stem is inserted, and the soft Wilson lines are 
modified as 
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(A.2) 



The offshellness is given by Ai = —pf/n-pi, and A2 = —v\l n 'Vi where pi (P2) is the n (n) 
collinear momentum of the collinear particles with the corresponding soft Wilson lines to be 
attached. Note that the insertion of Ai and A2 looks similar to the rapidity regulator, but 
it is the regulator for the IR divergence. Though similar in form, their sources are distinct. 
If we put the offshellness explicitly, Aj and 5i take different forms. In the soft Wilson line, 
A, can be obtained from the offshellness of a single collinear particle where the soft gluons 
are attached. On the other hand, the rapidity regulator <5, is obtained by the emission of the 
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n-collinear gluons from all the collinear or heavy particles in other directions. Therefore 
Si have complicated dependence on the offshellness of the other particles. Only in the 
back-to-back current, there exists a simple relation 8% = A 2 , 82 = Ai. 

In obtaining the dependence on the offshellness, we consider collinear particles or 
antiparticles from —00, or to 00, as considered in Ref. [20], and assign nonzero offshellness 
to the collinear particles. Here also arises the problem of gauge invariance. But it suffices 
to say that this is only an intermediate step to regulate IR divergences with the offshellness 
since the dependence of the offshellness is cancelled in the final results. 

Here we employ the dimensional regularization for the UV divergence, and the IR 
divergence appears as logarithms of Ai or A2. The virtual correction in Fig. 1 (a), and 
the real gluon emission in Fig. 1 (b) are given as 
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In calculating these matrix elements, the following plus distribution functions are used. 
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where Li 2 (x) is the dilogarithmic function. The one-loop result for the soft function in DIS 
is given by 
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Note that the soft function does not depend on A 2 , as it should be. In DIS, the final state 
is described by the final-state jet function which depends on A 2 , but A2 is not the IR 
cutoff. Instead it is related to the invariant jet mass of the final states. Therefore A 2 does 
not represent the IR divergence, and the soft function is independent of A 2 . 
For DY process, the soft function is defined as 
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The soft Wilson lines Y n and are defined in the same way as Y n and Yn except that n 
is replaced by n, and Ai by A2. 

The virtual correction and the real gluon emission are given as 
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In the final stage of computing M? DY , there are two possibilities in taking the limit 
Ai,A 2 — > 0. That is, the limit Ai — > can be approached first with A 2 fixed, and 
then the limit A 2 — > is taken. The limiting procedure can be reversed, however, the 
result is the same irrespective of the order of taking limits. 

The soft function in DY process obtained by adding the hermitian conjugate of Eq. (A.7) , 
and is given by 
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As can be seen again in Eqs. (A. 5) and (A.7), the soft functions contain IR divergences 
as well as mixed divergences. Therefore the soft functions themselves are not physical. 



B Collinear distribution functions with offshellness 

The collinear distribution functions can also be evaluated with the offshellness for the IR 
regulator. The poles in e are UV divergences. The naive collinear matrix elements are 
given as 

cx s Cf ^ rl / , (5\ , u 2 1, 9 5 , a 2 , S 7r 2 i 

O s Cf x / 1 , $ , J, II 2 1 , 2 5 7T 2 \ 

+ (T^ G " lnz + ln ^) - <TT^) + / 

M c = °^{l-z)(- + hx-^-2-\nz{l-z)). (B.l) 
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Here we put n ■ p = Q, and 5 is the rapidity regulator in the collinear Wilson line. If 5 is 
not used, IR poles in erR appear instead, and they also cancel. 
The soft zero-bin contributions are given as 



M, 
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2ir 
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Mj°) = 0. 



(B.2) 



Therefore the collinear contributions with the soft zero-bin subtractions are given as 
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(B.S) 



The collinear quark distribution function at one loop is given as 
f%{z, M) = 2Re (M a + M 6 ) + M c + 5(1 - z)(^ 1} + r£>) 



2tt 



(5 ( 1 - 2 )(4 + |- + - ln 

Ve z 2e e — p 



+ - In ^ + In 



(B.4) 
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2 + 2 6 J 
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1 + 



1 + z 2 
(1-z). 



In z - 2(1 - z) + (1 + z) ln(l - z) 



where is the counterterm and is the residue in the self-energy of the fermion £ n 



at one loop and they are given as 

4 47T £ ^ 47T V -p 2 / 



(B.S) 



The usoft zero-bin contribution differs only in the real gluon emission, which is given 



as 
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The corresponding collinear part with the usoft zero-bin subtraction is given as 

M Ms = M b - M, (0) 

o s Cf , / l 
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5(1- 



(B.7) 
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And the one-loop correction to the PDF is given as 
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The initial-state jet function can be written as 



2Re(- 
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To one loop, the kernel W are given as 





= 0, 
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(1-*)+ 

In this calculational scheme, the kernels are also IR finite. As in the case with the 
dimensional regularization, wjyjL = 0. However, wj^ Y 1S t ne same as the result with the 
dimensional regularization except the 7r 2 term. The source of the disparity can be seen 
from Eq. (A. 8). There is a term ln 2 (// 2 /p|) in Sdy, while there is ln 2 (p 2 / — p|) in K qq . 
Therefore the IR divergence cancels as expected, but there is a remnant of it 2 . It results 
in the difference of the term with ir 2 . If we take the limit p 2 ,p| — > 0, it corresponds to the 
dimension regularization limit for the IR divergence. Then the logarithms turn into poles 
in eiR, and there is no additional n 2 term involved. However, if we strictly keep the signs of 
the offshellness, this additional factor of ir 2 appears. If there were only single logarithms, 
this ambiguity does not occur. We have not been able to confirm whether the difference is 
due to the scheme dependence and further consideration is needed. 
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